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For intermediate Coulomb energy to Fermi energy ratios , spinless fermions in a random poten- 
tial form a new quantum phase which is nor a Fermi glass, neither a Wigner crystal. Studying small 
clusters, we show that this phase gives rise to an ordered flow of enhanced persistent currents, for 
disorder strength and ratios where a metallic phase has been recently observed in two dimensions. 
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An important parameter for a system of charged par- 
ticles is the Coulomb energy to Fermi energy ratio Ts . In 
a disordered two-dimensional system, the ground state 
is obvious in two limits. For large r^, the charges form 
a kind of pinned Wigner crystal, the Coulomb repulsion 
being dominant over the kinetic energy and the disor- 
der. For small Ts, the interaction becomes negligible and 
the ground state is a Fermi glass with localized one elec- 
tron states. There is no theory for intermediate , while 
many transport measurements following the pioneering 
works of Kravchenko et al and made with electron 
and hole gases give evidence of an intermediate metallic 
phase in two dimensions, observed for instance when 
6 < Ts < 9 for a hole gas in GaAs heterostructures. A 
simple model of spinless fermions with Coulomb repul- 
sion in small disordered 2d clusters exhibits a new ground 
state characterized by an ordered flow of enhanced per- 
sistent currents for those values of r^. In a given cluster, 
as we turn on the interaction, the Fermi ground state can 
be followed from Ts = up to a first level crossing. A sec- 
ond crossing occurs at a larger threshold after which the 
ground state can be followed to the limit 7's — s- oo. There 
is then an intermediate state between the two crossings. 
In small clusters, the location of the crossings depends on 
the considered potentials, but a study over the statistical 
ensemble of the currents supported by the ground state 



gives us two well defined values and 
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the system on a torus threaded by an Aharonov-Bohm 
flux, we denote respectively Ii and It the total longitu- 
dinal (direction enclosing the flux) and transverse parts 
of the driven current. One finds for their typical ampli- 
tudes |/(| ~ exp— (r^/rf) and /; « exp— (rg/r^) with 
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Below rf , the fiux gives rise to a glass of local 
currents and the sign of /; can be diamagnetic or para- 
magnetic, depending on the random potentials. Above 
rf , the transverse current is suppressed while an ordered 
flow of longitudinal currents persists up to , where 
charge crystallization occurs. The sign of /; can be para- 
magnetic or diamagnetic depending on the filling factor 
(as for the Wigner crystal), but does not depend on the 
random potentials (in contrast to the Fermi glass). One 
finds rf and r^ in agreement with the values delimiting 
the new metallic phase when 0.3 < kpl < 3, kp and / de- 



noting the Fermi wave vector and the elastic mean free 
path respectively. For kpl > 1, /; is strongly increased 
between rf and r^ ■ This suggests that the intermediate 
phase of our model is related to the new metal observed 
in two dimensions by transport measurement which we 
shortly review. 

In exceptionally clean GaAs/AlGaAs heterostructures, 
an insulator-metal transition (IMT) of a hole gas re- 
sults from an increase of the hole density induced 
by a gate. This occurs at « 35, in close agreement 
to « 37, where charge crystallization takes place ac- 
cording to Monte Carlo calculations , and makes highly 
plausible that the observed IMT comes from the quan- 
tum melting of a pinned Wigner crystal. The values of 
rg where an IMT has been previously seen in various 
systems (Si-Mosfet, Si-Ge, GaAS) are given in Ref. [||, 
corresponding to different degrees of disorder (measured 
by the elastic scattering time r). Those rs drop quickly 
from 35 to a constant value Tc ~ 8—10 when r be- 



comes smaller. This is again compatible with 
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given by Monte Carlo calculations g for a solid-fluid 
transition in presence of disorder. If the observed IMT 
are due to interactions, it might be expected that this 
metallic phase will cease to exist as the carrier density is 
further increased. This is indeed the case |^ for a hole 
gas in GaAs heterostructures at « 6 where an insu- 
lating state appears, characteristic of a Fermi glass with 
electron-electron interactions. 

In this work, we take advantage of exact diagonal- 
ization techniques for large sparse matrices (Lanczos 
method) where tiny changes of energy can be precisely 
studied. This restricts us to small clusters and low filling 
factors. Fortunately, the dependence on particle number 
has proved to be remarkably weak in many cases. In the 
clean limit, calculations [|j with 6-8 particles give the 
condensation of the electron gas into an incompressible 
quantum fluid when a magnetic field is applied. Pikus 
and Efros Q have obtained rY ~ 35 from 6x6 clusters 
with 6 particles, close to rY ~ 37 obtained by Tanatar 
and Ceperley for the thermodynamic limit. In the dis- 
ordered limit which we consider, there is another reason 
for expecting weak finite size effects. When the energy 
levels do not depend very much on the boundary condi- 
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tions, the periodic repetition of the same cluster cannot 
drastically differ from the thermodynamic limit obtained 
from an ensemble of different clusters. This usual local- 
ization criterion applies for insulators as the Fermi glass 
or the pinned Wigner crystal. Small cluster approxima- 
tions should then be sufficient for small and large . This 
explains why the critical factors which we will discuss 
are close to the thermodynamic limit given by the exper- 
iments. Finite size effects can be important only if one 
has a metal for intermediate r^. 

We consider a simple model of iV = 4 Coulomb in- 
teracting spinless fermions in a random potential defined 
on a square lattice with = 36 sites. The Hamiltonian 
reads: 

^ clc, +^«.n. + {/^|^. (1) 

<i,j> i i^j 

c| {ci) creates (destroys) an electron in the site i, t is the 
strength of the hopping terms between nearest neighbors 
(kinetic energy) and is the inter-particle distance for 
a 2d torus. The random potential of the site i with 
occupation number = cj^Cj is taken from a box distri- 
bution of width W . The interaction strength U yields a 
Coulomb energy to Fermi energy ratio = U / {2t^Tme) 
for a filling factor rie = N / L"^ . The disorder to hopping 
energy ratio W/t is chosen such that kpl takes values 
where the IMT has been observed . A Fermi golden 
rule approximation for r gives kpl ~ 192TTne{t/W)'^ . 
One has He — 1/9, W/t — 5,10,15 corresponding to 
kpl = 2.7, 0.67 and 0.3 respectively. 

The boundary conditions are always taken periodic in 
the transverse y-direction, and such that the system be- 
comes a torus enclosing an Aharonov-Bohm flux (j) in the 
longitudinal x-direction. Imposing (p = n/2 {(f) = tt corre- 
sponds to anti-periodic condition), one drives a persistent 
current of total longitudinal and transverse components 
given by 

and It = ^ /E respectively. The local current l\ flow- 
ing at the site i in the longitudinal direction is deflned 
by l\ ^ 2Im((«'o|c^^^^c,^_,J*o)) and by a correspond- 
ing expression for /*. The response is paramagnetic if 
J; > and diamagnetic if /; < 0. We begin by showing 
behaviors characteristic of a single cluster when varies. 

Fig. |l| corresponds to kfl < 1 (W/t — 15). Look- 
ing at the low energy part of the spectrum, one can see 
that, as we gradually turn on the interaction, classifica- 
tion of the levels remains invariant up to first avoided 
crossings, where a Landau theory of the Fermi glass is 
certainly no longer possible. Looking at the electronic 
density pi = (5'o|ni|\E'o) of the ground state |^o), we 
have checked that it is mainly maximum in the minima 



of the site potentials for the Fermi glass. After the sec- 
ond avoided crossing, pi is negligible except for four sites 
forming a lattice of charges as close as possible to the 
Wigner crystal triangular network in the imposed square 
lattice. The degeneracy of the crystal is removed by the 
disorder, the array being pinned in 4 sites of favorable 
energies. 




FIG. 1. Behavior of a single cluster for kpl < 1 as a func- 
tion of Vs. Top: Low energy spectrum (a 1.9rs term has been 
subtracted); Bottom left: jumps of 7 at the second crossing 
where the ground state (filled diamonds) and the first excited 
state (empty diamonds) are interchanged. Bottom right: /; 
(left scale, empty circles) and number of occupied sites 
(right scale, filled circles). 

For the same cluster, we have calculated C(r) = 
N~^^^PiPi-r and the parameter 7 = max^ C(r) — 
min^ C(r) used by Pikus and Efros Q for characteriz- 
ing the melting of the crystal. 7=1 for a crystal and 
for a liquid. Calculated for the ground state and the 
first excited state, ^{rg) allows us to identify the second 
crossing with the melting of the crystal. Moreover, one 
can see that the crystal becomes unstable in the interme- 
diate phase, while the ground state is related to the first 
excitation of the crystal (Fig. |^ bottom left). Around the 
crossings, the longitudinal current /; and the participa- 
tion ratio — N'^{J2i pf)^^ of the ground state (i.e. of 
the number of sites that it occupies) are enhanced (Fig. |l| 
bottom right). The general picture is somewhat reminis- 
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cent of strongly disordered chains ^ where level crossings 
associated to charge reorganizations of the ground state 
are accompanied by enhancements of the persistent cur- 
rents. Fig.|] is representative of the ensemble, with the 
restriction that the location of the crossings fluctuates 
from one sample to another as well as the sign (param- 
agnetic or diamagnetic) of /; below the first crossing, in 
contrast to Id. 




FIG. 2. Behavior of a single cluster for kpl > 1 as a 
function of rs- Left coordinates: currents /; (circle) and It 
(triangle). Right coordinates; (filled circle). 



Fig. H corresponds to kfl > 1 {W/t = 5). The previous 
level crossings are now almost suppressed by a stronger 
level repulsion and charge crystallization occurs more 
continuously. There is instead a broad enhancement of 
/; which, in contrast to Fig. ^, is not accompanied by a 
corresponding increase of ^s, which smoothly decreases 
from 20 of the 36 possible sites down to 4 when charge 
crystallization becomes perfect. A transition of the per- 
sistent current, from a disordered array of loops towards 
an ordered flow as increases has been noticed ||Tl[] by 
Berkovits and Avishai. To illustrate this phenomenon, 
the total transverse current It is shown in Fig. [H One 
can see that It is suppressed at rs ~ 5 while Ii continues 
to increase up to « 15. We have checked that a disor- 
dered array of loops persists up to « 5, followed by an 
ordered flow of enhanced longitudinal currents persisting 
up to w 15. The disordered array of loops gives rise 
to a diamagnetic or paramagnetic current depending 
on the microscopic disorder. The ordered flow gives rise 
to a paramagnetic However, Coulomb repulsions do 
not always yield a paramagnetic response. For instance, 
4x6 clusters with N = 6 become always diamagnetic 
at large r^. One can only conclude that the sign of the 
response in 2c? does not depend on the random potential 
when Vs is sufficient for suppressing It- In Id, Legett's 
theorem H states that the sign of /; depends on the par- 



ity of N only, for all disorder and interaction strength. 
The proof is based on the nature of "non symmetry dic- 
tated nodal surfaces " , which is trivial in Id, but which 
has a quite complicated topology in higher d. It is likely 
that such a theorem could be extended in 2d when the 
transverse flow is suppressed. 
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FIG. 3. Statistical study of an ensemble of clusters for 
W/t = 5, (circle) 10 (square) and 15 (triangle) as a func- 
tion of Ta. Top left: mean value < Ii >; Top right: fraction 
Cd of diamagnetic samples; Middle left: distribution of the 
logarithms of the paramagnetic current /;,p at = 1.7 and 
W/t = 15. Middle right: mean number of occupied sites 
by the ground state. Bottom left: longitudinal paramagnetic 
(empty symbols) and diamagnetic (filled symbols) currents. 
Bottom right: transverse currents. The straight lines are ex- 
ponential fits giving and rY shown in Fig. m 



We now present a statistical study of an ensemble of 
10^ clusters for W/t = 5, 10, 15. At the left top of Fig. |, 
one can see an increase of the mean // by about one order 
of magnitude when r-g « 7 for W/t — 5. We note that the 
persistent currents | [lO[ | measured in an ensemble of meso- 
scopic rings are typically higher by a similar amount than 
the theoretical prediction neglecting the interactions. At 
the right top of Fig. |^, the fraction of diamagnetic clus- 
ters is given as function of rs , showing that the enhance- 
ment of the mean is partially related to the suppression of 
the diamagnetic currents. This suppression is faster for 
weak disorders. The mean number of sites occupied 
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by the ground state is given at the middle right of Fig. 
~|, showing a neghgible increase when W/t = 15 at low r^. 
and a regular decay otherwise. The paramagnetic Ii^p and 
diamagnetic Ii^d longitudinal currents, and have log- 
normal distributions for all values of rs when W/t > 5. 
The stronger is the disorder, the better is the log-normal 
shape of the distribution (see middle left of Fig. ^) . The 
average of the logarithms give the typical values shown 
in the bottom part of Fig. ^. On the left, the longitudi- 
nal currents // are given, the diamagnetic responses Ii^d 
(filled symbols) being separated from the paramagnetic 
responses Ii^p (empty symbols), while the transverse cur- 
rents It are given at the right side. The log-averages 
exponentially decay as Ii^d oc \It\ oc exp— (r^/rf) and 
Ii^p oc exp— (r^/r^) when is large enough. The vari- 
ances of log|/t| and log// increase as rg/rf and r^/r^ 
above rf and respectively. The values of rf and 
extracted from the exponential fits (straight lines of 



Fig. 1^) are given in Fig. 
diagram is proposed. 



where a sketch of the phase 
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FIG. 4. Proposed phase diagram for 2d spinless fermions in 
a random potential. rY (filled circle) and rf (empty circle) 
obtained from Fig. ^ bottom. 

Fig. H and Fig. ^ show that a simple model of spinless 
fermions with Coulomb repulsion in a random potential 
can account for the critical carrier densities and disor- 
der strengths where the IMT occurs. The comparison 
between the curve ^'^(t) given in Ref. |^ (summarizing 
the factors where the IMT has been observed) and 
the curve rY {kpl) of Fig. ^ (characterizing the suppres- 
sion of /;) is very striking. The value Vs — 6 where the 
reentry has been observed in Ref. |^ is also compatible 
with the curve {kpl) characterizing the suppression of 
It- We have not indicated in the proposed phase dia- 
gram the difference between kpl > 1 (where /; has a 
strong enhancement which may be the signature of a new 
metal at the thermodynamic limit) and kpl < 1 (where 
/; persists up to without noticeable enhancement). 
A study of the size dependence at a fixed filling factor 



will be necessary for studying a possible IMT driven by 
an increase of kpl at intermediary r^. This is unfortu- 
nately out of reach of exact diagonalization techniques. 
Another striking difference is that and Ii convey sim- 
ilar information when fc^^ < 1 while the increase of /; 
is accompanied by a decrease of when kpl > 1. This 
suggests that transport for intermediary results more 
from a collective motion of charges than from a delocal- 
ization of individual charges. The spin degrees of free- 
dom are not included in our model, the orbital part of 
the wave function being totally anti-symmetrized. This 
restriction is quite important for short range screened 
interactions, but is certainly less severe for long range 
interactions and low densities. However, there are many 
experimental evidences [|2| that spin effects play a role. 
A more complex phase diagram is possible when spins 
are included. But even for 2d spinless fermions, these 
small cluster studies allow us to conclude that there is 
a new quantum phase, clearly separated from the Fermi 
glass and from the Wigner crystal, identified by a plastic 
flow of currents without charge crystallization, which is 
likely to give a new metal in the thermodynamic limit 
when kpl > 1. 
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